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JOSEPHSON  JUNCTION  TRIANGULAR  PRISM 
QUBITS  COUPLED  TO  A  RESONANT  LC  BUS 

Qubits  and  gates  for  a  holonomic  quantum  computer 


S.  P.  Yukon 

Air  Force  Research  Laboratory,  Electromagnetics  Technology  Division,  Hanscom  AFB, 

MA  01731,  USA 

Abstract:  We  investigate  the  properties  of  Josephson  junction  triangular  prism  qubits  coupled  by 
mutual  inductance  to  a  resonant  LC  bus.  We  show  how  the  symmetries  of  the  qubit  potential  may  be 
used  to  implement  the  Duan  Cirac  Zoller  scheme  for  holonomic  quantum  computation. 

Keywords:  Josephson  junction  qubit.  Geometric  gates,  Holonomic  quantum  computation 


1.  INTRODUCTION 

A  geometric  approach  for  achieving  quantum  computation  using  non-Abelian 
holonomic  gate  operations  has  been  recently  proposed  by  Pachos  et  al.  [1],  An 
implementation  of  (non-Abelian)  holonomic  quantum  computation  using  a  trapped  ion 
quantum  computer  [2]  has  subsequently  been  proposed  by  L.  M.  Duan,  J.  I.  Cirac.  and 
P.  Zoller  (DCZ)  [3].  We  propose  here  an  alternative  physical  implementation  of 
holonomic  quantum  computation  based  on  a  previously  described  quantum  computer 
employing  Josephson  junction  (JJ)  triangular  prism  qubits  coupled  to  a  resonant  LC  bus 
[4].  For  the  ion  trap  system  the  excitation  of  different  transitions  can  be  distinguished  by 
polarization  or  by  frequency.  For  qubits  based  on  the  symmetric  potential  well  of  a  JJ 
triangular  prism,  symmetric  and  antisymmetric  excitations  of  prism’s  two  JJ  triangles 
using  external  magnetic  flux  take  the  place  of  laser  pulse  polarization.  Because  the  qubit 
potential  well  and  qubit’s  interaction  with  external  flux  are  defined  by  Bohm  Aharonov 
vector  potential  integrals  around  each  cell  circuit  (yielding  an  Abelian  holonomy  or  Berry 
pha.se),  the  basis  for  both  qubit  and  gate  properties  can  be  said  to  be  of  geometric  origin. 

The  DCZ  qubits  are  atomic  ions  that  have  the  desirable  feamre  that  all  of  the  qubits 
possess  identical  properties  when  they  are  isolated  from  local  environmental  perturbations. 
The  JJ  qubits  that  we  propose  to  employ  have  inherent  fabrication  errors  and  no  two  qubits 
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will  be  exactly  identical.  To  overcome  this  variability,  the  local  microwave  rf  sources  that 
provide  the  gate  pulses  at  each  qubit  need  to  have  their  frequencies  adjusted  to  match  the 
transition  frequencies  needed  for  that  particular  qubit.  A  desirable  aspect  of  holonomic 
quantum  computing  for  such  nonidentical  qubits  is  that  once  the  local  rf  source  frequencies 
have  been  adjusted,  it  is  possible  to  cany  out  a  set  of  global  (or  multi-local)  logic  operations 
as  the  holonomic  gates  depend  only  on  relative  phases  and  amplitudes  of  the  rf  gate  sources 
and  not  on  their  absolute  amplitudes.  Because  the  gate  rotation  angle  for  a  holonomic  gate 
depends  on  the  solid  angle  subtended  by  the  closed  path  in  parameter  space,  holonomic  gates 
should  be  relatively  immune  to  small  random  fluctuations  in  external  magnetic  flux  since  the 
integral  of  the  fluctuations  they  induce  on  a  closed  path  should  average  to  zero. 

1.1  Josephson  junction  triangular  prism  qubits 

In  Fig.  1  we  have  sketched  the  three  alternative  versions  of  the  flattened  JJ  prism  qubit 
proposed  in  Ref.  [4].  Of  these  the  four-  and  five-cell  versions  have  similar  properties 
except  that  the  five-cell  version  has  an  extra  constraint  condition  (on  the  sum  of  the  fluxes 
through  the  five  cells)  that  gives  rise  to  unwanted  combinations  of  second-order 
excitations.  Of  the  two,  the  simpler  four-cell  version  is  thus  preferable. 

The  three-  and  four-cell  versions  may  be  distinguished  by  the  circulating  currents  in 
the  outer  cells  when  the  system  is  at  a  potential  well  minimum.  The  expressions  for  the 
effective  one-dimensional  potential  for  both  the  three-  and  four-cell  versions  may  be 
written  as  Vrix)  —  — 2cos(y)  -t-  rcos(2y)  (shown  in  Fig.  3  for  r  =  0.675).  The  expression 
for  the  phase  of  a  (dotted)  junction  at  a  potential  minimum  in  one  of  the  triangular  cells 
may  be  written  as  ^*  =  ±x*  +/(4',)  where  ±y*  are  the  two  minima  for  Vrix),  and/(fl>,) 
represents  a  function  of  fl>,s  (where  is  the  external  flux  in  the  ith  cell).  The  circulating 
current  in  the  left-  and  right-hand  JJ  triangular  cells  at  the  two  minima  may  be  represented 
as /mean  ±  A/ (left-hand  cell)  and /mean  +  A/ (right-hand  cell).  For  the  three-cell  circuit 
(with  fl),.  =  0  defined  below)  /mean  =  0.415,  A/  =  0.415  at  r  =  1,  and  /mean  =  0.65, 
A/ =  0.32  at  r=0.65.  For  the  four-cell  circuit  (with  <1)^  =  2fl>^  defined  below) 
/mean  =  0.0,  A/  =  0.87  at  /■  =  1,  and  /mean  =  0.0,  A/  =  0.65  at  r  =  0.65.  The  four-cell 
circuit  yields  a  larger  A/  for  gradiometer  coupling  to  a  SQUID  detector  and  is  therefore  the 
circuit  that  we  choose  for  the  qubit. 

The  computational  basis  states  for  the  qubit  are  taken  to  be  the  ground  and  first  excited 
state  of  the  potential  Vr(x)-  Because  the  ground  and  first  excited  states  wavefunctions  of  a 
symmetric  potential  are  purely  symmetric  and  antisymmetric,  respectively,  there  is  an  equal 
probability  of  the  current  being  in  the  right-  or  left-hand  JJ  triangle.  To  perform  a 
measurement,  a  Hadamard  gate  or  a  7r/2  T-rotation  on  a  qubit  in  one  of  the  basis  states  will 
lead  to  a  superposed  wavefiinction  that  is  centered  around  the  left  or  right  potential  well 


Figure  1.  Sketch  of  three  versions  of  the  flattened  triangular  JJ  prism. 
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minima.  This  in  turn  will  lead  to  one  of  the  previously  described  current  patterns  that  can  then 
be  measured  by  a  SQUID  detector.  A  measurement  perfomied  without  the  Hadamard  gate 
would  be  able  to  measure  symmetric  or  antisymmetric  combinations  of  the  qubit  basis  states. 

While  we  have  used  the  symmetric  properties  of  the  potential  V^ix)  to  define  the  qubit, 
we  have  used  the  results  of  Mooij  et  al.  with  regard  to  the  low  inductance  properties  of  the 
qubit  [5-7]  and  the  necessity  for  suppressing  transitions  to  other  charge  states  [5, 6].  We  will 
not  discuss  the  aspects  of  coupling  the  qubit  to  detectors  and  the  environment  here,  but  the 
results  of  Refs  [5]  and  [6]  have  addressed  these  issues  and  are  applicable  to  the  qubit 
proposed  here. 

The  DCZ  Hamiltonian  employs  three  ground  or  metastable  states  plus  a  common 
excited  state  for  each  ion  qubit.  The  Hamiltonian  in  the  rotating  frame  for  the yth  ion  with 
rf  sources  for  each  needed  transition  turned  on,  is  given  by 

Hj  =  h^lo\e)j{0\  +  bn^\e)j{l\  +  hili,\e}j{b\  +  h.c.  (1) 

where  11,  are  the  Rabi  frequencies  for  the  transition  from  state  |()  to  state  \e).  Transcription 
of  the  DCZ  holonomic  quantum  computation  scheme  to  Josephson  junction  qubits 
requires;  (1)  qubits  with  three  degenerate  states  (in  the  rotating  frame)  that  couple  to  a 
common  excited  state;  (2)  a  way  of  exciting  each  of  the  three  states  without  exciting  the 
others  (and  being  able  to  control  the  amplitude  and  relative  phase  of  the  separate  excitation 
modulation  envelopes);  and  (3)  a  bus  having  harmonic  oscillator  energy  levels  that  can 
couple  to  all  of  the  qubits.  In  the  next  two  sections  we  demonstrate  how  this  may  be 
accomplished  using  four-cell  JJ  prism  qubits  coupled  to  a  resonant  LC  bus. 

The  Hamiltonian  (in  units  of  2Ej)  for  the  four-cell  qubit  in  Fig.  2,  derived  in  Ref.  [4], 
may  be  written  as 

=  -  I 2i7//  cos cos (x  +  Y^+P cosilx) cos( </>,  +  |  (2) 


Figure  2.  Flattened  Josephson  junction  prism  qubit  with  and  S<(),  control  lines. 
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where  we  have  defined  6  =  0]  —  4>,/2  and  A"  =  •Ai  ~  <Ao)/2  with  =  (4>,ii  + 

<At/2)/2,  along  with  the  diagonalizing  +  combinations;  lAi  =  (fi  —  ^2)/2.  4^2  =  (fi  “  &)/2, 
^1  =  (^1  +  f2)/2,  02  —  (li  +  f2)/2.  as  well  as  (A.,.,,  =  (<A„i  +  <t>a2)!2.  with  the  definitions 
<A<,,.r.«  =  2ir  ,„/<!>(),  Ej  =  /f4’()/2'n' and  $0  =  h/2e.  We  have  used  the  zero  inductance 
limit  Bohm-Aharonov  phase  constraint  equations  {{]  —  ^2  “  ^2  =  <Aai.  ~  ^2  ~  ^2= 
‘A«2-  ~  +  ^2  =  <A</2>  ^2  “  =  ‘A</i }  for  the  gauge  invariant  and  ^Junction  phases.  For 

junctions  ^1,  ^2-  and  ^2-  the  junction  capacitances  and  critical  currents  have  been  taken 
equal  to  C  and  Ic,  respectively.  For  junctions  ^3  and  ^3,  the  junction  capacitances  and 
critical  currents  have  been  taken  equal  to  pC  and  pic  respectively,  allowing  the  mass  term 

to  be  written  as  —  (2Ej)M  =  (1  +  2p)  (2Ej)/{4(Ec/2)),  with  Ec  =  e  ^l{2C).  The 
T]fj  term  comes  from  the  Hartree  replacement  cos(0)  ^  (O«|cos(0)|Oh)  =  rj/y.  For  the 
ground  state  10^),  p/y  1.  In  order  to  introduce  the  terms  of  the  Hamiltonian  that  can  be 
utilized  for  qubit  gates,  we  further  divide  (A^y,  <A«.  and  <A.,  into  constant  and  small  excursion 
constituents  as  <A^  =  (A^  +  <A„  =  <A^  +  5<A„,  and  <A,  =  <AI^  +  S<A,.  We  take  <Ar  =  0  ('-o. 

no  constant  or  time-varying  antisymmetric  component  of  the  flux  threading  the  middle 
two  cells). 

For  this  Hamiltonian,  the  time-independent  Schrodinger  equation  is  equivalent  to  the 
Whittaker  Hill  equation,  for  which  exact  solutions  have  been  found  [SJ.  The  lowest  even 
and  odd  solutions  are  denoted  by  gCo(A'),  SC2(x)^  and  gs^ix)  and  are  defined  in  terms 

of  continued  fractions.  If  we  identify  the  computational  basis  states  (;t'|l)  and  (;t'|0)  with 
gC(}ix/2)  and  gs\{x/2),  the  assignment  of  gC2(A'/2)  and  gs-^ix!'^)  to  the  DCZ  auxiliary 
states  and  (;t'|/2)  will  allow  all  of  the  DCZ  one  and  two  qubit  holonomic  gates  to  be 
carried  out.  The  two  basis  states  and  auxiliary  states  are  shown  in  the  inset  of  Fig.  3  along 


Figure  3.  Energy  levels  and  qubit  transitions  (inset)  for  V,.);^')  (in  units  of  lEj)  with  r  =  0.675. 
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with  the  complete  set  of  bound  state  levels  for  the  qubit  potential  of  Eq.  (2)  written  as 
V,(.X)  =  ‘-'os(  X  +  -  rcosdx)  with  r  =  -pcos(</)"  +  4^l)/{r]„  cos((^||/2))  =  0.675. 

Since  we  have  chosen  to  use  the  excited  states  gcjlxd)  and  gST.lxJ'^)  as  the  DCZ 
auxiliary  .states,  their  energy  eigenvalues  (and  that  of  gc^jx!^)^  which  is  a  possible 
transition  from  g^2(A'/2))  must  lie  appreciably  below  the  potential  well  saddle  that 
separates  unit  cells  in  the  two-dimensional  periodic  potential  (in  order  to  prevent  the 
excitation  of  different  charge  states)  [5, 61.  There  is  a  range  of  values  of  and  rj/y  for 

which  this  is  possible.  As  an  example  for  the  choice  r  =  0.675,  =  100,  the  saddle  lies 

0.74  X  (2£/)  above  the  potential  minimum,  while  the  gc4,{xl'^)  uud  gSi(xl2)  states  lie 
respectively  at  0.18  x  (2£/)  and  0.12  x  (lEj)  above  the  minimum. 

A  second  requirement  for  the  potential  is  that  the  transition  frequencies  (U2i  = 
(£2  ~  E\)/h  and  0)43  =  (£4  —  £3)//?  be  separated  by  at  least  10  x  Xl/„  where  O,,  is  the 
Rabi  frequency  for  the  gs\(xl2) gc^lx/'^)  transition,  in  order  not  to  excite  the 
gsy{x/2)  ^  gC4{x/2)  transition. 

To  model  the  inclusion  of  rf  forcing  and  coupling  to  other  qubits,  we  expand  the 
potential  in  a  power  series  around  =  0,  0,  and  4>^  #  0,  keeping  Scj),  terms  up  to 

second  order.  Physically,  the  and  magnetic  flux  can  be  supplied  by  a  uniform 
external  field.  The  excitation  can  be  supplied  by  a  vertical  line  bisecting  the  qubit 
while  the  Sc/),,  excitation  can  be  supplied  by  two  symmetric  horizontal  control  lines 
designed  to  null  or  minimize  concomitant  excitation  of  Sc/)^  and  5</),.  as  shown  in  Fig.  2 
(the  wanted  Sc/),,  excitations  will  come  from  the  bus).  The  rf  excitations  ,,  are  tuned  to 
the  transition  frequency  between  \e)  and  |0),  1 1 ),  or  |S).  Thus  the  full  expression  for  S(/)„  for 
the  S  e  transition  is  thus  given  by 

5</)<;-*>  =  (3) 

where  84>^l’\t)  is  the  adiabatically  varying  modulation  of  the  rf  pulse  envelope.  The  time 
varying  phase  may  be  introduced  as  a  frequency  modulation  of  the  carrier  at  the 
transition  frequency  {E,.  —  Ei,)/h  with  the  FM  pulse  frequency  given  by  = 

(£,,  -  Ei,)/h  +  d<p,,(t)/dt. 

The  Hamiltonian  for  the y'th  qubit  in  the  rotating  frame  with  all  rf  sources  on,  may  thus 
be  expressed  as  the  DCZ  Hamiltonian  of  Eq.  (1)  with  the  Rabi  frequencies  for  the  three 
needed  transitions  given  by 


0(1  =  r7„(0|sin(;^')|e')cos(<^^/2)S(/);'’* 

(4a) 

and 

ill  =  p(l|cos(2;k')l£>)sin((/)‘’  -|-  </)‘))S</)*,'’ 

(4b) 

=  '9//(/^|sin(;^)k)co.s((/)“/2)5</);f’ 

(4c) 

1.2 

Coupling  the  resonant  LC  bus  to  the  qubits 

To  enable  two-qubit  gates  to  be  cairied  out,  each  of  the  qubits  is  coupled  to  a  resonant 
LC  bus  by  mutual  inductance  coupling  between  a  loop  in  the  bus  and  the  two  inner  cells  of 
the  qubit  as  shown  in  Fig.  5.  The  bus  is  based  on  the  trapped  ion  bus  scheme  of  Refs. 
[3,9, 10]  and  is  similar  to  the  resonant  LC  bus  coupling  of  qubits  in  Refs.  [11-13].  The 
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operators  for  flux  and  charge  in  an  LC  oscillator  with  resonant  frequency  w  =  1  js/LC,  can 
be  written  in  terms  of  creation  and  annihilation  operators  {a'*',  a}  as 


<1>  = 


(a  +  a  ), 


(5) 


If  there  are  N  qubits  coupled  to  the  bus,  and  a  fraction is  available  for  coupling  to  the 
qubits,  the  amount  of  flux  coupled  to  each  qubit  will  be/y^h/M  The  second-order  terms  in 
the  Hamiltonian  expansion  will  then  yield  interaction  terms  for  even-even  (or  odd-odd) 
transitions  fV|<7^’|(a -f  a ' )  and  even-odd  transitions  iVQoj|j(a  4- a  ’ )  and  i^tTj),(a  +  a^), 
where 


and 


(tJi^  =  |c)^(;u,|  +  h.c.  Hermitian  conjugate 


(6) 


n;, 


(0|sin(A')k)277f, 


^hLui/l 


pcos(<^"  -I-  <^“)(l|cos(2A')|<?)277f, 


yfhLwp. 


<a|sin(A')k)277f^ 


s/hLo}/2 


(7a) 

(7b) 

(7c) 


Adopting  the  Sprensen  and  Mplmer  [5, 6]  bichromatic  bus  scheme  used  by  DCZ,  which 
requires  only  virtual  bus  excitations,  the  eflective  Hamiltonian  for  a  pair  {j,  fc}  of  qubits 
coupled  by  the  resonant  LC  bus  may  be  written  as 


2Ej 


(8) 


This  expressions  differs  from  those  in  DCZ  and  Refs.  [5, 6]  due  to  the  fact  that  we  have 
taken  into  account  the  slightly  different  parameters  of  each  qubit,  and  have  allowed  for  a 
third  detuning  frequency  offset  So  in  addition  to  the  DCZ  frequency  offset  5. 

The  total  capacitance  C  will  consist  of  the  capacitance  between  the  wires  forming  the 
bus  and  any  additional  capacitance  added.  The  bus  frequency  m  =  is  chosen  to 

avoid  direct  excitation  of  transitions  in  the  qubit  and  to  optimize  the  number  of  qubits  in 
the  system.  To  remain  within  the  quasi-static  limit,  the  maximum  extent  of  the  bus  is 
~0.1A,v,  where  A^.  is  the  substrate  wavelength.  Assuming  a  bus  frequency  of  3  GHz,  and  a 
substrate  refractive  index  n,  =  3.0  yields  A,/ 10  =  0.33  mm.  Assuming  (lEj)/ 
=  200  GHz,  with  r  =  0.675,  100  the  expressions  (4c)  and  (7c)  with  = 

0.001  X  2-77  yield  n/,/27r=  0.152  GHz  and  lT^/27r=  0.627/\/AMHz.  The  Rabi 
frequencies  for  two-qubit  gates  given  in  Eq.  (8)  may  be  increased  from  this  value  by 
taking  the  offset  frequencies  5  and  So  to  be  small  fractions  of  HJ,,  and  by  choosing  n  such 
that  the  rf  frequency  Wrf  =  (£2  —  E\)lh  +n-  o)  is  distant  from  any  allowed  transition 
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frequency  so  that  the  value  of  84>^,  could  be  increased.  An  example  would  be  to  choose  n 
such  that  ftjri  ^  (£7  —  E\)/h. 

We  note  that  by  setting  (f)^  =  0,  the  middle  cells  in  Fig.  1  may  be  eliminated  and 
their  junctions  combined.  The  bus  for  these  two-cell  versions  should  then  be  moved  to 
the  right  (or  left)  cell.  From  Eqs  (4)  and  (7)  we  see  that  for  Fig.  1(b),  all  one-qubit 
teiTTis  would  survive  but  only  the  even  two-qubit  term  (7b)  would  survive  for  =  0. 
This  would  still  allow  a  universal  set  of  gates  to  be  formed  but  without  the  extra  degree 
of  control  from  the  even-odd  transition  terms  (7a)  and  (7c).  The  two-cell  version  of 
Fig.  1(a)  would  have  similar  lo.sses.  For  Fig.  1(c),  all  even  transition  terms  vanish  [14]. 
For  one-qubit  gates,  an  even  term  could  be  fa.shioned  for  |1)  ->  \e)  by  using  second- 
order  odd  transitions  since  the  energy  denominator  for  the  intermediate  state  transition 
|1)  ^  |0)  is  small. 


1.3  Non-Abelian  holonomic  gates 

The  basic  set  of  DCZ  single  qubit  gates  f/W  =  g'hdj' ^ 

qubit  gate  a^g  sufficient  for  all  of  the  operations  needed  to  carry  out 

holonomic  quantum  computations.  The  DCZ  Berry  phases  <^|,  ipy,  and  4>n,  are  given  by 
their  implicit  dependence  on  6  and  cp  as: 

= —^^sin{d)dd  dtp  (9) 

with 

n,,  =  0,  ii,  =  -n'  sin/^^^c'>,  life  =  n"  cos 


4>y  —  osin(6) do  dtp 


with 

Oo  =  fisin(0)cos(<p),  fl|  =  n'sin(6)sin(<p),  fl/,  =  fl"cos(0) 
with  (1  =  il'  —  n"  and 


(10) 


<t>\\ 


— -()sin(6)  d(p 


(11) 


with 


Du  =  0, 


=  tan(? 


The  minimal  set  of  gates  for  quantum  computation  (two  single-qubit  gates  and  one  two- 
qubit  gate)  must  form  a  noncommuting  set,  which  implies  that  their  geometric  phases  must 
be  generated  by  non-Abelian  holonomies  [1].  The  computational  basis  during  an  adiabatic 
evolution  cycle  rotates  into  “dark”  states  that  maintain  zero  eigenvalues.  An  example  of 
dark  states  for  the  gate  Hamiltonian  would  be  |Dyi)  =  cos(ff)(cos(^)|0)y -1- 

sin(^)|l)y)  —  sin(<p)|^>)y  and  \Dy2)  =  cos(9)|0)y  —  sin(^)ll)y.  Since  the  gate  Hamiltonians 
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have  zero  eigenvalues  for  these  states,  there  is  no  dynamical  phase  contribution  from  the 
gates  to  the  computational  basis  states  [1]. 

Using  expressions  for  the  geometric  phase  in  [1],  the  DCZ  set  can  be  completed  with 
the  single-qubit  X  rotation  and  two-qubit  and  rotations  where  and  cTy  are 

replaced  by  {|a;S)<a'/3'l -1- |a/3)<a'/3'l}  and  {-(lai3Xa')3'| -h  t|a'i8')(a/3|}.  respectively, 
and  |a/3)  (and  la'jS'))  is  one  of  the  two-qubit  states  {|0,0a),  |0;U),  11/0*),  |1,U)}.  The 
expression  for  the  X  rotation  is  given  by  where 


=  c)sm(6)dd  d(f> 


(12) 


with 


fl()  =  flsin(0)cos(^),  fl|  =  n'sin((l)sin(9)e'""^",  fl/,  =  fl"  cos(  ) 
For  Yap.a'p’  rotation  the  two-qubit  Hamiltonian  can  be  written  as 


+  in' ||n*|cr«',TS 

rotation  is  then  given  by 


The  Y, 


(13) 


with 


with 


(fiy  ^  sin{0)dd  d(p  and 


ilnp  =  Osin(ff  )cos(9), 


(14) 


^a'p'  —  fX  sin(0)sin((p),  11/,  =  n"cos(6) 

The  Xap.a'^  rotation  can  be  similarly  derived  following  the  pattern  of  the  single-qubit  X 
and  Y  rotations. 

1.4  Fault  tolerance  of  Ui  and  Un  and  Uy  and  holononiic  gates 

Assuming  that  it  will  not  be  possible  to  adjust  the  three  microwave  sources  at  each 
qubit  to  obtain  equal  basic  Rabi  frequencies  il  for  the  three  needed  transitions,  we  have 
looked  at  the  effect  of  unequal  basic  Rabi  frequencies  {fl,  fl',  fl")  on  the  value  of  the 
resulting  gate  rotation  angle.  We  define  an  “input”  gate  rotation  angle  </>j„  as  the  Berry 
angle  given  in  Eqs  (9),  (10),  or  (1 1 )  for  a  perfectly  adjusted  system  with  a  three-step  pulse 
sequence:  {0=O^7r/2,  (p=0},  [9='tt/2,  ^=0^^iinai},  {(1=  77/2^0,  (p  = 
‘Piinail-  The  pulse  envelope  used  is  0(0  =  ^(tanh(A(r -|- /o))  —  tar>h(A(t  —  fo)))/2, 
<p(f)  =  <p(tanh(Af  -h  l)/2,  with  A  =  1/2,  fo  =  20,  0  =  7r/2,  il  =  Ha  =  27r.  The  gate 
rotations  presented  below  are  calculated  in  the  restricted  Hilbert  space 
7^0  =  {|0),  |1),  \e),  |0)}.  In  Fig.  4(a)  we  have  plotted  the  resulting  gate  rotation  angle 
for  the  (/,  gate  for  {H,  H',  O"}  =  Ha{0,  1.0,  1.0),  HaIO,  1.0,  1.2),  and  Ha{0,  1.0, 
1.4).  Since  all  of  the  gate  rotation  angles  lie  on  a  straight  line  for  a  given  set  of  basic  Rabi 
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Figure  4.  (a)  Output  vs.  input  rotation  angle  for  basic  Rabi  frequency  error  sets  0,  1,  2;  (b)  Output 

vs.  input  rotation  angle  for  basic  Rabi  frequency  error  .sets  0,  1,2,  3. 


frequencies,  measurement  of  a  single  pair  of  input  and  output  rotation  angles  will  then 
determine  the  needed  input  angle  for  any  desired  output  rotation  angle.  The  behavior  of  all 
gates  like  U\  \  that  have  flo  =  0  will  be  similar. 

In  Fig.  4(b)  we  have  plotted  the  output  rotation  angle  for  the  f/Km,,  gate  for  {fl,  O', 
O"}  =Oa{1.0,  1.0,  1.0},Oa{1.0,  1.22,  0.83),  Oa{  1.00,  1.44,  0.66),  and  Oa{1.00,  0.66, 
1.34).  For  this  gate  all  of  the  output  rotation  angles  for  a  given  set  of  basic  Rabi 
frequencies  lie  on  a  45°  rotated  sine  curve.  Measurement  of  a  single  pair  of  input  and 
output  rotation  angles  can  then  determine  the  needed  input  angle  for  any  desired  output 
rotation  angle.  The  behavior  of  all  gates  like  Ux,  Uy,  and  that  have  adiabatic 

envelopes  similar  to  Eq.  (10)  will  have  behavior  similar  to  that  of  V .  Unlike  the  [/|  i 
and  U\  gates,  there  are  also  two  points  at  ±  tt/I  where  the  output  rotation  angle  perfectly 
matches  the  input  rotation  angle.  Thus  any  of  the  gates  Ux.  Uy,  Ux^i^j^y,  or  at 

rotation  angle  ±  tt/I  will  be  fault  tolerant  with  respect  to  microwave  pulse  amplitude 
errors.  Investigation  of  the  fault  tolerance  properties  of  gates  when  the  Hilbert  space  is  not 
restricted  to  Tfo  requires  optimizing  the  potential,  the  resonant  pulses,  and  bus  parameters 
to  minimize  leakage  out  of  Tfo,  and  will  be  addressed  in  a  future  paper. 

1.5  Connecting  buses  into  a  network 

It  is  also  possible  to  connect  groups  of  buses  in  an  open  branching  network  as  first 
described  in  Ref.  [4].  To  transcribe  this  feature  to  the  holonomic  qubits  and  gates 
described  above  requires  distinguishing  a  particular  qubit  (or  qubits)  that  is  to  be  assigned 
the  sole  task  of  transferring  gate  excitations  from  one  bus  to  a  second  and  vice  versa.  This 
is  shown  schematically  in  Fig.  5.  Since  the  f/rij  ,,  is  fault  tolerant  at  rotation  angle  +7r/2, 
we  shall  employ  it  to  connect  buses  and  define  it  as 


Ucy,A±)  =  (lOAXOAl  +  |0/U.><o,i,|  +  n.uXiAi  ±  UAXXUI}  (i5) 

To  carry  out  a  two-qubit  t/cr,j  (  + )  gate  between  a  control  qubit  j  on  bus  1  and  a  target  qubit 
k  on  bus  2  via  the  common  transfer  qubit  t  (assumed  to  be  in  its  ground  state),  it  is  necessary 


158 


S.  P.  Yukon 


Figure  5.  Schematic  of  two  coupled  buses  and  transfer  (center)  qubit. 


to  first  carry  out  a  t/cr,,(— )  gate  with  qubit  j  the  control  qubit  and  qubit  t  the  target  qubit.  A 
second  Ucy,j,(—)  gate  with  qubit  t  the  control  and  qubit  k  the  target  and  a  final  UcyjA~)  gate 
with  qubit  j  the  control  qubit  and  qubit  t  the  target  qubit  achieve  the  de.sired  result  with 
control  qubit  t  returned  to  its  original  ground  state.  While  they  are  not  perfectly  fault 
tolerant,  three  control-not  gates,  between  ( j,  t),  (f,  k),  and  ( j,  t)  would  also  have  the  overall 
effect  of  a  control-not  between  j  and  k  with  t  reset  to  its  original  ground  state. 

To  carry  out  a  control-y  over  a  number  of  buses,  it  is  necessary  to  first  carry  out  a 
sequence  of  Ucyj(-)  starting  from  the  control  qubit,  through  all  of  the  intervening  transfer 
qubits  to  the  target  qubit,  followed  by  a  final  sequence  of  t/cn.,(~)  from  the  original  qubit 
to  the  last  transfer  qubit  in  order  to  reset  the  transfer  qubits  to  their  original  ground  states. 
The  sequence  involving  two  transfer  qubits  would  thus  be  (/ct,/— )  between  (y, /i), 
(ti,  Ui,  k)  followed  by  resetting  transfer  qubits  with  Ucy,j{—)  applied  to  (ti,  tj)  and 
( j,  f I ).  The  maneuvers  described  in  the  preceding  paragraph  may  be  chained  across  as 
many  coupled  buses  as  required  (given  an  environment  where  perfect  coherence  is 
possible).  This  allows  open  branching  networks  of  coupled  buses  to  be  formed,  enlarging 
the  possible  number  of  interacting  qubits  in  a  single  quantum  computer. 


2.  CONCLUSIONS 

We  have  shown  that  it  is  possible  to  design  a  system  comprised  of  Joseph.son  junction 
prism  qubits  coupled  inductively  to  a  resonant  LC  bus  that  is  capable  of  carrying  out 
holonomic  quantum  computations.  The  effective  Hamiltonian  is  zero  throughout  an  entire 
computation,  so  there  is  no  dynamical  evolution  of  the  system.  In  the  restricted  Hilbert 
space  Kti,  one-  and  two-qubit  X  and  Y  rotation  gates  are  shown  to  be  tolerant  of  pulse 
amplitude  errors  for  rotation  angles  of  +  7r/2.  The  behavior  of  one-  and  two-qubit  gates  at 
arbitrary  rotation  angles,  when  pulse  amplitude  errors  are  present,  can  be  predicted  if  one 
measurement  can  be  made  to  determine  an  appropriate  input/output  curve. 
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